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We develop a fully fledged theory of quantum dynamical patterns of behavior that are nonlocally 
induced. To this end we generalize the standard Laplacian-based framework of the Schrodinger 
picture quantum evolution to that employing nonlocal (pseudodifferential) operators. Special at- 
tention is paid to the Salpeter (here, m > 0) quasirelativistic equation and the evolution of various 
wave packets, in particular to their radial expansion in 3D. Foldy's synthesis of "covariant particle 
equations" is extended to encompass free Maxwell theory, which however is devoid of any "particle" 
content. Links with the photon wave mechanics are explored. 



I. INTRODUCTION 
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] The standard unitary quantum dynamics and the Schrodinger semigroup-driven random motion, are examples 

1^ ' of dual evolution scenarios that may be mapped among each other by means of a suitable analytic continuation in 
Ph I time procedure. This is an offspring 0-0 of Euclidean quantum field theory methods, albeit reduced to the purely 
\^ I quantum mechanical level. Both evolutions are generated by means of a common local Hamiltonian operator. Our 
departure point for subsequent analysis is an observation that a complete spectral resolution of the corresponding 
Hamiltonian actually determines a classical space-time homogeneous diffusion- type Markov process in R" . 

Within the general theory of so called infnitcly divisible probability laws (see below) the familiar Laplacian (Wiener 



■7* ', noise or Brownian motion generator) is known to be one isolated member of a surprisingly rich family of non- 



Gaussian Levy noise generators. All of them stem from the fundamental Levy-Khintchine formula, and typically 
[ refer to probability distributions of spatial jumps and the resultant jump-type Markov processes. That needs to be 
^ . contrasted with the traditional diffusion imagery (Wiener noise and process) associated with the Laplacian, Q . 
O^' The emergent Levy generators are manifestly nonlocal (pseudo-differential) operators and, while being additively 
perturbed by a suitable external potential, give rise (via a canonical quantization procedure described subsequently) 
I ' to Levy- Schrodinger semigroups. The dual image of such semigroups comprises unitary dynamics scenarios which can 
^ , be viewed as signatures of a nonlocal quantum behavior. As we discuss in below this dynamical nonlocality extends 
00 ' to the very concept of photons, within so-called photon wave mechanics, 

C ■ I Quite apart from "Euclidean" vs "real" time connotations, the considered dual dynamical systems refer to real 
\l i time labels and clocks. Both of them drive probability density functions (pdfs) which, in the course of time evolution 
" I ' {t S i?+) either maintain or develop asymptotic heavy-tails and typically have a finite number of moments in existence. 

, This needs to be contrasted with the Gaussian standards of thinking (all pdf moments in existence, rapid decay at 

' infinities etc.) that pervade the Laplacian-based quantum theory. 

. The major goal of the present paper is to set on solid grounds the quantization programme that completely avoids a 
. . ' reference to classical mechanics (of massive particles), normally viewed as a conceptual support for the specific choice 
^ I of the Hamiltonian operator within the traditional Schrodinger wave mechanics. Here we consider a standard form of 
the Hamiltonian (minus Laplacian plus a perturbing potential) as an exception rather than a universally valid feature 
of an admissible quantum theory. 
, The latter is introduced by means of the most primitive quantization ansatz, whose core lies in choosing the Hilbert 
space L^(i?") as an arena for our investigations. From the start we have the Fourier transformation realized as a 
unitary operation in this space and a canonical quantization input as a straightforward consequence. 

The above mentioned Levy-Khintchine formula, while being tailored to our purposes, derives through a Fourier 
transform of a symmetric probability density function. A variety of symmetric probability laws for random noise is 
classified by means of a characteristic function which is an exponent T](p) of the (27r)"/^-multiplied Fourier transform 
of that pdf: / d"xp(x) exp{zLipx) = exp[r](ji)]. 

The most straightforward, naive canonical quantization step p ^ p = — iV (hereby limited to the symmetric pdf 
case) introduces the notion of random transport that is driven by Lcvy-Schrodinger semigroups cxp[tri{p)]. Their dual 
partners actually are the unitary evolution operators cx'p[itri{p)] of interest and are a subject of further discussion. 

In the present paper we develop a complete theory that starts from diffusion type processes and their semigroup 
reconstruction, followed by the duality mapping. Next we pass to general jump-type stochastic processes that stem 
from the Levy-Khintchine formula and pass to the dual (e.g. quantum) dynamics, with a number of specific physically 
motivated examples, like the Salpeter Hamiltonian and its m = Cauchy version. The major departure point for the 
present analysis is Ref. Q , but earlier references of relevance can be traced back therein, see also Q . 

Our main objective is to deduce and next expose a fully nonlocal pattern of dynamical behavior appropriate for the 
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description of general quantum dynamics. We believe that it is nowadays necessary to take under scrutiny customary 
ways of thinking about quantum phenomena and open new conceptual avenues, based on the inherent nonlocality of 
the dynamics generator, c.f. [l^l for an earlier attempt. 



II. DIFFUSION-TYPE PROCESSES VS QUANTUM MOTION 
A. Ground state condition 



Taking as obvious the standard wisdom about the Schrodinger picture quantum dynamics, we merely recall (here, 
without specifying the space dimension) that the equation 

i' = (1) 

involves a Hamiltonian H that is a self-adjoint operator in a suitable Hilbert space domain; we consider time- 
independent potential V = V{x). Since we are interested in non-negative generators of motion, for later convenience 
we impose an additive renormalization of a bounded from below Hamiltonian (to be absorbed in the explicit functional 
expression for V), such that H >Q and is its lowest eigenvalue. 

Further standard notations are reproduced for the record: for any normalized ip^ Pi^jt) = I^P(2;,i) stands for a 
pdf, V = {h/2mi)[{\/ip/ip) — {Vip* is a current velocity field, while j = vp a probability current: dtp = — Vj. 
While keeping the polar (Madelung) decomposition of in mind , we consider v in the gradient form: v ~ ^^^s. 

Wc note that by introducing u(x,t) = (?i/2m) Vlnp and next the drift field b{x,t) = u + w, we can rewrite the 
continuity equation in the Fokker-Planck form dtp = D Ap — \7{bp), with D = playing the role of the diffusion 
coefficient. Here, in general b = b{x, t) is time-dependent, unless we pass to stationary states. 

1/2 

The ground state condition for H corresponding to the lowest eigenvalue 0, i. e. Hil) = with ji/'l — pj , directly 
involves the negative of the familiar de Broglie-Bohm "quantum potential" 

V = 2nW' ^ . (2) 

P* 

Here D = h/2m is implicit, although not a must if out of the quantum context (like in the case of Brownian motion). 

1 /2 

Having , we have completely determined V = which in turn allows for a complete reconstruction of the 

(Schrodinger) semigroup cxp{—tH /2mD), [O, A traditional well established procedure actually amounts to the 
reverse: given the (Schrodinger) semigroup with a priori chosen (admissible) potential V{x), one infers pl^'^[V] = (hj"^ 
as the lowest strictly positive eigenfunction of the semigroup generator (i.e. of —H/2mD). 



ihdtip 
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A + V 



B. Schrodinger semigroup vs Fokker-Planck dynamics 

We keep the notation D for a diffusion coefficient. The (Schrodinger) semigroup induces the related generalized 
diffusion equation (by setting V ~ we get the standard heat equation) 



exp{-tH/2mD)'^o = * =^ ^i* = 



2mD 



^ = — * (3) 

2mD ^ ' 



This Schrodinger-type equation (albeit without an imaginary unit i) employs H as & self-adjoint relative of the 
Fokker-Planck operator, 

Lfp = DA~V{b-) =^ dtp = Lfp p (4) 

and thence indirectly determines the time evolution of a pdf p[x, t) for t > 0, the initial data being tacitly presumed. 

To execute a transformation between the Fokker-Planck and the Schrodinger semigroup dynamics, let us assume 
that the pertinent FP pdf p{x,t) has a strictly positive asymptotics p*{x). We make an ansatz 

^i^,t) = Jp±. (5) 



3 



Additionally, wc identify b = DVlnp* as a timc-indcpcndent forward drift, and introduce (typically time-dependent) 
osmotic u = DV In p and current v = b — u = (1/to)Vs velocity fields. 

The connection between the Fokker-Planck and semigroup dynamics is readily established, provided a compatibility 
condition Eq. (2) is extended to the drift field b = DV In 



Vix) = 2mD 



2 



(6) 



It is clear that ^{x,t) solves the generalized diffusion equation (3) if the identity (5) is valid. Given V{x), we can in 
principle solve the Riccati-type equation (6) with respect to b{x). Proceeding in reverse, with a pre-defined forward 
drift b{x) of a diffusion-type process, we have determined an admissible form of the Schrodinger semigroup potential 
V{x). 

We point out that in the theory of diffusion processes, an impact of external conservative forces is typically encoded 
in the drift function b ^ —WU where U is a. Newtonian potential. We emphasize a substantial difference of functional 
forms for V and U (except for the notorious harmonic oscillator case). 



The preceding discussion r nay be summarized by a symbolic identity, appropriate for operators acting in their 

-1/2 J- 1/2 _ - jj /»N 

-p* LfpP* =- — -H. (7) 



n r nay 

domains. Namely, we have P, 

1 



2mD 

where HpJ = 0. On the other hand, if we consider H as given a priori, then the formula Lpp ~ —pj' 
associates a corresponding Fokker-Planck generator Lpp with the Hamiltonian H . 
It is worthwhile to mention the equivalence of Eq. ([7]) and the transformation recipe: 

dtp = D [py'A(p-^/^ . ) - p:"\/\pl'^)] p (8) 

which directly follows from Eq. (3), once we set p = vj/ p^^, see also [l3| . 

This formula is particularly interesting, since it effectively transforms the free noise (Wiener) generator into its 

" . 1/2 

additive perturbation, e.g. iJ, provided that p* is known. This observation will prove useful in the Levy flights 
context. In fact Eq. has the structural form (up to the replacement of A by any Levy noise generator) of the 
general transport equation for p(a;, t) that encodes a response of non-Gaussian noise to external potentials, see below. 

C. More on Lfp H transformations 

Before passing to non-Gaussian jump-type processes, few comments should be made to bridge notational and 
conceptual incongruences between the well rooted methodology of the statistical mechanics approach to random 
motion, [l[ and math-oriented relevant references @, [l^ . 

In the theory of Markovian diffusion processes, physicists pay a particular attention to the Fokker-Planck equation 
and thence Lpp, while mathematical practitioners advocate an operator Lpp = DA + 6V that is Hermitian adjoint 
to Lpp. Actually, L*pp is an infinitesimal diffusion generator that is directly involved in the uniqueness proof for the 

Markov process in question. We note that Eq. ([7]) implies —pl^^ L*ppp^ = H/2mD 

If we have in hands a transition probability density p{y, s, x,t), s < t, a, necessary condition for the existence of a 
unique Markovian dynamics is that p solves the so-called backward (first Kolmogorov) diffusion equation 

dsP^ -DAyp-bWyP= -L*pp{y,s)p, (9) 

where b — b{y, s). The very same transition probability density solves also the forward (second Kolmogorov) equation 

dtp = DA,p-V{bp)^Lpp{x,t), (10) 

where b — b{x, t). The latter (forward) equation, once rewritten with respect to a pdf p[x, t) ~ J p{y, sq, x, t)p(y, so)dy, 
with So < t, is named the Fokker-Planck equation by physics practitioners. 

We indicate that p{y, 0, a;, t) Ax needs to be interpreted as a probability that a particle trajectory started from y 
at time will reach a Aa; vicinity of a point x at a time t > 0, e.g. J p{y, s, x, t)dx = 1 for all s < t and all y E R. In 
the mathematical literature, p{y, 0, x, t) = p{t, y, x) = pt {x\y) would typically appear. 
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More formal description of the above comments, that is close to a standard mathematical lore, looks as follows, 
[T^ . Given a Markov transition function Pt{x\y) of a random process on R. The generator of the process L reads 

1 r+oo 

(Lf){y) = lim-[/ pMy)f{^)dx - f{y)] (11) 

If a transition function is stochastically continuous (a priori presumed, see e.g. [3): then the corresponding semigroup 
Tt is defined by 

/ + 00 
Pt{x\y)f{x)dx (12) 
-oo 

and is strongly continuous, so that its generator L is densely defined in a suitable domain. 

In such a case we can also define an adjoint semigroup T^* acting on the space of (probability) densities L^{R, dx), 

/ + 00 
ptix\y)p{y)dy (13) 

We denote its generator by L* . In case of diffusion- type processes we obviously have the generic outcome L* = Lpp 
and L = L*pp. 

All that may be checked by inspection in case of the familiar Ornstein-Uhlenbeck (OU) process, We scale away 
physically relevant constants taking L — b{y\/y + {1/2) Ay as a process generator, where b{y) = —y. The OU transition 
density reads 

p{y, 0, x, t) = pt{x\y) = [^(1 - e-2*)]-i/2 exp 

Here L* = (1/2)A^ — '^x\b{x) ■] and h{x) = ~x. This process is stationary (time- homogeneous), hence we can safely 
replace thy [t — s) > 0, thus passing to p{y, s, x, t) and the related forward and backward (Kolmogorov) equations. 

Remark : One may explicitly resort to dimensional units, with x being replaced by the velocity label v, @. Then 
the transition probability density of the OU process takes the form 

with s < t and has an asymptotic {t — )• oo) invariant density p{v) = (27r/3D)~^/^ • exp{—v^/2f3D). The drift of the 
process reads b{v) ~ —f3v and p solves the Fokker-Planck (second Kolmogorov) equation dtp = DA^p — V„(6p). 
Here /3 is a friction coefficient and D = kBT/m/3 reficcts Einstein's fiuctuation-dissipation relationship; ks is the 
Boltzmann constant, T stands for an equilibrium temperature of the thermostat. 



{x 



(l-e- 



(14) 



D. Dynamical duality: "Euclidean" mapping exemplified 

In the light of the above discussion, compare e.g. Eqs. (1) and (3), it appears quite persuasive to execute (at least 
formally) the Wick rotation in the complex time plane it — > t > 0; ?i — > 2mD, under the restriction t e i?+ : 

eyi-p{-iHt/h)i;o = =^ exp(-tif/2mD)*o = *t (15) 

that maps between diffusion-type and quantum mechanical patterns of dynamical behavior, [IHli- In the absence of 
external potentials (free case) we have H = —2mD^A. 

For clarity of discussion, it is instructive to invoke explicit examples. We pass to one spatial dimension and rescale 
(or completely scale away) a diffusion coefficient. Given a spectral solution for H = —A + V > in L^{R), the integral 
kernel of ex.p{—tH) reads {t -t- it gives rise to the kernel of cxp{—itH)) 



k{y,x,t) = k{x,y,t) ^Y.'^M-^jt) *j(y)*K^)- 

3 



(16) 
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Remember that we assume eo = and the sum may be replaced by an integral in case of a continuous spectrum. 
Then one needs to employ complex-valued generalized eigenfunctions, like e.g. ^j{x) — > ^p{x) = (27r)~^ exp(ip2;). 
Indeed, if we set V{x) — identically, a familiar heat kernel is readily obtained : 



k{y,x,t) = [cxp(iA)](y,x) 



1 

2^ 



cxp(— p^t) cxp(ip(y — x) dp 



(17) 



(47ri)-i/2 exph(2/-a;)V4i], 

in accordance with / exp(— cr^p^) exp(— ipa;)(ix = (tt/ct)^/^ exp(— p^/4(T^), where cr > 0. We note that the kernel of 
[exp(t_DA)] appears after changing the timc-scalc in ([TT]), t —J- Dt. A formal identification D = 1/2 gives the kernel 
that is often met in the mathematical literature and corresponds to (1/2) A instead of A. 

Consider H = (l/2)(— A + x^ — 1) (e.g. the rescaled and (—1) renormalized harmonic oscillator Hamiltonian). The 
integral kernel of cxp(— ti?) is given by a rescaled form of the classic Mchlcr formula, (Tsl. [l6j: 



1 °° \ 

Hy,x,t) = [exp(-tiJ)(y,x) = -j= c^p[-{x^ + y^)/2] ^ ^Hn{y)Hn{x) e^p{-nt) 



(18) 



(^[l-exp(-2t)])-i/2 



exp 



1/2 2n (a^-e 'y? 

-^i- -y)- (i_,-2.) 



where e„ — n, $„(a;) = [4"(n!)^7r]~^/'* exp(— a;^/2) _ff„(a;) is the L'^{R) normalized Hermite (eigen) function, while 
Hn{x) is the n-th Hermite polynomial Hn{x) = (— l)"(expx^) exp{~x^). 

The normalization condition J fc(y, x, i) exp[(y^ — x'^)/2]dy = 1 actually defines a transition probability density of 
the Ornstein-Uhlenbeck process (see e.g. Eq. ((TO)) ) 



p{y, 0, X, t) = p{y, X, t) = k{y, x, t) pV^ [x) / pV^ {y) 

with p*(x) = TT"-^/^ exp(— x^). 

A more familiar form of the Mehler kernel reads (note the presence of exp(t/2) factor) 



(x^ + y^) coshi — 2xy 



2 sinh t 



(27rsinht)i/2 

By formally executing t ^ it one arrives at the free Schrodinger propagator 

K{y,x,t) = [exp(i<A)](2/, x) = (27r)~^/^ / exp(— zp^t) exp(ip(y — x) dp = 



(19) 



(20) 



(21) 



(47rit)"^/^ exp[+i{y ~ xf/U] 
and likewise, at that of (—1 renormalized) harmonic oscillator propagator 

. (x^ + y"^) cost — 2xy 



K{y, X, t) = /^^P('^/2) 
(27ri smi)^/^ 



2 sint 



(22) 



Learn a standard Euclidean (field) theory lesson concerning multi-time correlation functions. In the exemplary 
harmonic oscillator case, t > t' > 0; t ^ it results in: 



E[X{t')X{t)] = / p^{x')x' p{x',t',x,t)xdxdx' ^ (1/2) cxp[-{t-t')] 



(23) 



W{t',t) = {i;o,qH{t)qH{t')4'o) = (1/2) expH(t - t')] 

where qnit) stands for the position operator in the Hciscnberg picture, P. IT^. 

The major message of our discussion is that we have encountered the dual dynamical patterns of behavior that 
follow equally realistic clocks. The Euclidean mapping (Wick rotation) is merely a mathematical artifice transforming 
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one dynamical model into another, [1, H, with an evolution following a common time scale. Obviously, once on 
the quantum level, we may extend the validity of the formalism from t G i?^ to all times t €z R. 

Remark 1: In relation to the "Euclidean" label, we may invoke the statistical physics lore of the 50-ies and 
60-ties, by passing to an integral kernel of the density operator, that is parameterized by equilibrium values of the 
temperature. To this end one should set e.g. t = hu/kBT for a harmonic oscillator with a proper frequency uj and 
remember about evaluating the normalization factor 1/Zt, where Zt stands for a partition function of the system. 

Remark 2: In the confined (due to a suitable choice of an external potential) regime, the ground state cj), Hcj) = 0, 
of a nonlocal self-adjoint and non- negative Hamiltonian-type operator H, induces a pdf |0p = p*, which actually is 
an asymptotic target of an associated jump-type stochastic process. The initial knowledge of the ground state of H 
and the wellposedness of the affiliated random motion problem both ensure the existence, and allow for a complete 
reconstruction of the nonlocal Lcvy-Schrodinger quantum dynamics and its dual semigroup partner, see for example 



III. LEVY (JUMP-TYPE) PROCESSES AND NONLOCAL RANDOM DYNAMICS 

A. Levy-Khintchine formula 

Let us consider a family of infinitely divisible probability laws and affiliated stochastic Markov processes X{t) — ?> 
p{x,t) that are characterized by the celebrated Levy-Khintchine (LK) formula (for any spatial dimension n > 1) 

{exp[ipX{t)]) := / exp{ipx)p{x,t) dx = exp[—tF{p)] (24) 



R 



where 



/ + 00 
[expiipy) - 1 - -^]u{dy) (25) 

and the integral in Eq. ([25]) is interpreted in terms of its Cauchy principal value (that is often made explicit by 
denoting /fj„^Q instead of the "plain" integral J); ^{dy) stands for so-called Levy measure such that 

" u{dy) < oo (26) 



(1+2/^) 



A functional form of the Levy measure is uniquely singled out by specifying the Levy-Khintchine exponent F(p) 
and in reverse. We shall confine further attention to the family of symmetric stable laws associated with F^ij)) = 

where /i € (0,2) and Fm{p) = ^ p^ + m? — m,m > 0. Here, to eliminate inessential dimensional parameters, we 
have chosen suitable units. We note that we have in fact enco untered a re- scaled version of the classical relativistic 
Hamiltonian, which is better known in the dimensional form ^y rn^c^ + c^p'^ — m(? ^ where c is the velocity of light. 

Its m ^ limit is well defined and equals c\p\ which directly refers to the Cauchy (symmetric stable [i — X) 
probability law. We point out as well that Eq. represents a reduction of the complete LK formula to the integral 
term which determines jump-type probability laws and processes. We have disregarded a Gaussian contribution of 
the form F(p) = p^ /2, identifying the Wiener noise and process, whose properties underly the reasoning of Section II. 

To see the problem from a broader perspective, we recall that a characteristic function of a random variable X 
completely determines a probability distribution of that variable. If this distribution admits a density we can write 
i?[cxp(ipX)] = Jjj, p{x) C'xjp{ipx)dx. Infinitely divisible probability laws are classified by the general Levy-Khintchine 
formula 

E[eM^pX)] = exp[-i^(p)] = exp{zap - (aV2)p2 + f [exp{ipy) - 1 - -^]v{dy)} (27) 

J-oo 1 + y 

where v{dy) stands for the Levy measure. 

To elucidate the meaning of the deterministic term, we invoke the Cauchy probability density function in the form 
p{x) = {(j/t:)[{x — a)^ 4- cr^]~^. Its characteristic function reads i?[exp(ipX)] = exp(iap — <7\p\) and the term — (t|p| 
can be verified to come out from the (scaled by u) integral contribution in Eq. (p7)) . 

A direct consequence of an infinite divisibility of the considered probability laws is that the "free" noise fully 
determines a corresponding Markov process through 

E[eM^pXt)] = cxv[-tF{p)] . (28) 
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Accordingly, any composite "free" Markov process X{t) = Xt can be interpreted (decomposed) as follows: Xt = 
at + aBt + Jt + Mt where Bt stands for the free Brownian motion (Wiener process) , Jt is a Poisson process while Mt 
is a general jump- type process (more technically, martingale with jumps). 

In particular, by disregarding the deterministic and jump-type contributions in the above, we are left with the 
Wiener noise X = aB. For a Gaussian pdf p{x) ~ {27ra^)^^/^ exp{—x^ /2a^) we directly evaluate E[exp{ipX)] = 
exp(— CT^p^/2). It is enough to set = 2Dt to recover the induced Wiener process with E[exp{ipX{t)] = exp{—tDp^) 
and F{p) = Dp'^. 



B. Canonical quantization 



At this point we employ a substitution procedure, that is in fact a canonical quantization step, up to the explicit 
presence of h: 

p^p^ -iV =^ F{p) F{p) . (29) 

Wc point out that, in view the standard Fourier representation in use, a casual quantum mechanical operator notion 
{xf){x) = xf{x) is implicit. No covariant position operator, like that of the Newton- Wigner type, jl8|, is hcrc-by 
addressed. 

The domain of F{p) is tacitly placed in i^(i?"), whose element is in fact any considered p{x). That, irrespective of 
its L^{R^) normalization. A direct consequence of the L^{R'^) arena choice for the action of all considered operators 
is that the Fourier transform is a unitary operation. 

This is sufficient for the validity of standard position-momentum type uncertainty relations, since they automati cally 
follow from the properties of the direct and inverse Fourier transforms. This should be compared with Refs. [l^, l20l |. 
where Hardy's theorem concerning an interplay between the localization of a function and that of its Fourier transform 
has been analyzed (a folk transcript reads: a function and its Fourier transform cannot both be both very small). 
Here we also refer to Ref. [2l| for a discussion of uncertainty relation in relativistic quantum mechanics and Rcf. 
[2^ [2^ for a complementary view upon uncertainty relations in random motion. See also (2^ . [25j . 

If F{p) is given in the Levy-Khintchine form ([27)1 . then exp[— tF(j5)] is a contractive semigroup operator. In 
particular this pertains to exTp{tDA) where A denotes the spatial Laplacian and the Hamiltonian {l/2mD)H = 
F{p) = Dp^ = —DA, has been previously employed in the discussions of Section II. Note that we can get rid of the 
constant D by rcscaling the time parameter. 

From now on we shall restrict considerations to symmetric probability distributions of pure random jumps. The 
(infinitesimal) Levy measures v{dy) are odd with respect to the spatial inversion, e.g. i'{—dy) = —v{dy). Therefore, 
the regular form of the Levy-Khintchinc formula (|27p . while reduced to the integral expression (j25p . becomes further 
simplified (Cauchy principal value of the integral is implicit, see e.g. (sl. H^. [28l| ) : 

p+CO 

F{p) = - / [exp{ipy) - l]v{dy) (30) 



which in view of the canonical quantization recipe (|29p defines the action of the semigroup generator —F{p) on 
functions in its domain according to: 

- F{p)f{x) = +l [fix + y)- f{x)]dv{y) . (31) 

We emphasize that a generically singular behavior of the Levy measure in the vicinity of zero needs the (counter)term 
containing —f(x) for consistency reasons. There is no clean way to eliminate this contribution from relevant integral 
formulas. This prohibitive statement should be kept in memory while getting through a number of physics-oriented 
papers on related topics, as examples sec e.g. formulas (2.11) through (2.26) in [s^ or (14), (15) in Ref. [s^, c.f. the 
next subsection. 

As mentioned before, in the family of infinitely divisible probability laws our attention is focused on symmetric 
stable laws associated with F^{p) = \p\^^ where ^ g (0,2) and Fm{p) = y/jP+rri? — m with m > 0. Since many 
of our arguments will not rely on the specific spatial dimension we shall reproduce n-dimensional versions of the 
corresponding Levy measures. In particular, the ^ G (0, 2)-stable Hamiltonian operator if^ = F^,(p) = |A|^/^, 
induced by F^{p) = acts upon functions in its domain as follows [43ll43|: 



- |Ar'/V(.) = / [/(. + y) - m] Mdy) - ZZKl f 



W2|r(-f)| J \x-y 



d^'y (32) 
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where x € i?" and v^{dy) stands for a (self-defining) Levy measure ^ All above integrations are understood 

in the sense of the Cauchy principal value. 

Here, -|A|^/2 is a nonlocal (pseudo-differential, fractional) generalization of the ordinary Laplacian A. It is 
associated with symmetric stable probability laws and related jump- type stochastic processes (here called Levy flights), 
3 El, |33]-[i3l- We note that the coefHcient in Eq. ([32]) has been chosen to secure /^„[1 — cos{py)]dv^{y) = \p\^, 
43| . compare e.g. also Eq. (pO]) . 

The relativistic (named also quasi-, semi-, pscudo-relativistic or Salpeter) Hamiltonian operator H^^ — F„^{p) = 
[\/—A + — m], induced by Fm{p) = \/p'^ + rn? — m,m > 0, acts in its domain according to 26|: 

^ [T^aT^- m]/(.x) = I ifix + y) fix)] .M) = 2 (£) J [/(y) - /(.)] ^'^ ■ (33) 

It is instructive to notice that the to — !■ limit of the relativistic semigroup generator exists and coincides with /i = 1 
(Cauchy) stable generator — |A|. In the above K^{z) is the modified Bessel function of the third kind of order i^. 



C. Levy flights: transition pdfs and transport equations. 

"Free" semigroup kernels coincide with transition probability densities of the pertinent jump- type processes. As 
a complement to definitions of the previous two subsections, let us introduce calculational tools that are based on 
Fourier transformation techniques. An advantage of a functional analytic lore is that contractive semigroup operators, 
their generators and the pertinent integral kernels can be directly deduced from the Lcvy-Khintchine formula, see 
also 0. 

We extract from semigroup generators the "free" Hamiltonians, defined up to dimensional constants which can be 
easily recovered if needed. They have the form Hq = F{p), where p = — iV stands for the momentum operator (here 
we put H = 2mD = 1). Let f{x) be function in the domain of F{p) and f{p) = (27r)~"/^/ f{x)exp{—ipx)(Px its 
Fourier transform. We define 

/ + 00 
exp(-ti^(fc)) cxp(ifcx)/(fc)d"fc = [exp(-tF(p))/(p)]^(a;) (34) 
OO 

where the superscript V denotes the inverse Fourier transform. Let us set 

fct = (2^)-"/2[exp(-t^^(p)]^. (35) 

Then the action of exp(— tifg) can be given in terms of a convolution (i.e. by means of an integral kernel kt = 
k{x -y,t) ^ k{y,0,x,t)): 

eM-tHo)f = [exp{-tF{p))f{p)]'' ^ f * h (36) 

where (/ * g){x) ;= J gix — z)f{z)dz is a convolution of functions / and g. 

In view of Eq. ([55)1 we have in hands Fourier integral redefinitions for all pertinent semigroup kernels. We note 
that those kernels play the role of transition probability densities of the underlying Markovian jump-type processes. 
In particular, k^(x — y,t — s) has the form 



k^{x, t) = (27r)-" / exp(-ipx - t\pn d^p (37) 

while kmix ^ y,t — s) reads 



k^{x, t) = (27r)"" J exp[-ipx - t(yjm'^ +p'^ - to)] d> . (38) 
For completeness one should recall a Fourier form of the heat kernel k{x — y,t — s), compare e.g. Eq. P7)) . 

k{x, t) = [exp(<A)] (x) = (27r)-" / exp(-ipa; - tp^) d''p , (39) 



All the above (generalized) heat kernel formulas stem from a spectral resolution of the corresponding (rescaled) 
Hamiltonian operator, compare e.g. our discussion of Eqs. (|16p . ([TT]) and 
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One should be aware that closed analytic outcomes of Fourier integrals are scarce in the present framework. We shall 
reproduce explicit formulas for integral kernels of the Cauchy and quasi-relativistic semigroups, see e.g. [isl [26l. [27j. 
Namely, upon reintroducing suitable physical constants {h is kept equal 1, see however [30|), the integral kernel of the 
semigroup operator exp[— t(v^c?A + m^c? — mc^)] reads: 

fc„(x,t) = 2 (£y"''''^'c("+3)/2[i eMnicH)]{x' + c2i2)-("+i)/4 /f(„+i)/2(mcV^2 + cH^) ■ (40) 

In view of t~^km{y, t)dy — )■ Vm{dy)as t after putting c = 1, we retrieve the previous formula for the Levy measure 
as employed in Eq. (|33| . 

The mass m — > limit of fc,„(x,i) does exist and equals to the Cauchy kernel ko{x,t): 

ko{x, t) = r (41) 

whose more popular form (in view of computational simplicity) employs c = 1. 

In view of the semigroup dynamics rule [exp(— iiJ)/](x) = f{x, t), we readily get transport equations for pdfs driven 
by "free" noise. The equations 

dtp{x,t) = -\Ar/'pix,t) = -\Vrp = j[p{x + y,t)-p{x,t)]v^{dy) (42) 

where /i € (0, 2) and 



dtp{x, t) = - [ V- A + to2 - m]p{x, t) = [p{x + y,t) - p{x,t)]v,nidy) (43) 



are so-called master equations for the pertinent jump-type processes and, in the present setting, replace the standard 
Fokker-Planck type equation, appropriate for a diffusive transport. We recall that the Cauchy principal value need 
to be attributed to the integrals involved and one must not disregard the ^ p{x) counter term in those formulas. 

As an exemplary formula we reproduce the master equation for the ID Cauchy transition probability density 
ki{x — y,t — s) = p{y, s, x, t) with s < t: 

dtp{y,s,x,t) ^ [ [p{y,s,x + z,t) - p{y,s,x,t)]diyi{z) = [ ^^^''^'^'^^ — x, t) 
J J 7r(x - zY 



D. Continuity equation dtp = — Vj for symmetric stable processes. 



1. Inversion of S/ . 



In case of the Brownian motion (Wiener process) the continuity equation dtp ~ ^Vj is merely another form of the 
heat equation dtp = DAp, where a definition of j = vp along with v = — DV Inp does the job. Things become less 
obvious if the free transport equation is to take the Levy form appropriate for a symmetric stable process. 

For concreteness, instead of dtp ~ — 7|V|^/9 = —"f\A\'^^'^p = — 7Vj, with p e (0,2) and 7 being the jump intensity 
parameter, let us discuss the Cauchy case p = 1 in some detail. We note that m?c^ + c^p"^ — mc? in the mass m = 
limit reduces to c\p\. Accordingly we may pass to the nonlocal equation dtp = — c|V|p, c.f. Eq. (|37p which has a 
classic Cauchy solution 

p{x,t) = ki{x,t) = — / cxp(-ifca; - |fc|A) dfc = -— (45) 

2n J_^ ttV + x^ 

with A = \{t) — b + ct, 6 > 0. We note also that 

1 r+°° \2 _ ^2 

\V\p=— \k\eM-^kx-\k\X)dk= (46) 

27r7_oo 7r(A^+a;^)^ 

Accordingly dtp = —c\\/\p ~ — Vj, where an analytic expression for the probability current j{x,t) is given in terms 
of an indefinite integral (a symbolic inversion V^^ of V): 



(47) 
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Here v{x,t) = j/p = cx/X and we formally have cV~^|V|p = j- In view oi X = b + ct, ior sufficiently large times 

V ~ x/t holds true. We can also write v = — (c/2A)Vlnp. This should be compared with the Brownian outcome 

V = — DV Inp '-^ +x/t (to this end set e.g. p ^ cxp{—x^/2Dt)). 

The inversion formula V^^Vj — >■ j reported in Eq. (|47p may be interpreted in terms of the Fourier representation. 
Namely, the Fourier transform j{p) of j{x,t) reads 

If.. c 

j{k,t) — / cxp{~ ikx)j{x,t)dx — —isgnik) cxp(— A|fc|) . (48) 

V Stt j v27r 

We define the action of V^^ through (cannot be extended beyond ID): 

V""'^ exp(— ifcx) = exp(— i/cx) . (49) 

{-ik) 

Accordingly, we have from Eq. (|46p (|fc|/A: — sgnik), next employ k — —p) 

c\/~^\\/\p{x,t) ~ —— / sgn{k)exp{—ikx — \k\X) dk ^ / i{p,t) exp{ipx)dp = j{x,t) (50) 

One must be aware that rather restrictive assumptions are necessary for a function g to belong to the domain of 
V^^. In fact we can use only when g is "sufficiently nice". Such phrase is sometimes used in the mathematical 
literature in the analogous context of properly defining |A|~''/^ as a legitimate operator, [l^- As shown above, a 
function g = V/ where / is sufficiently fast decreasing at infinities, is an example of a "sufficiently nice" function. 



2. Inversion of \A\'^^'^ . 

One of our motivations is to analyze the origin and meaning of formulas of the type (— A)''/^~^V$ appearing in 
the literature devoted to so-called fractional quantum mechanics, also in relation to the probability current concept, 
(isl - lsoj . There- in, under the restriction 1 < p < 2, hence excluding the previously considered case of /.t = 1 and 
generally p € (O7 !)• Another motivation is related with the Landau-Peierls notion of the photon wave function, the 
definition of which employs (— A)^^/"', see e.g. (tI. [Sll - lssj . We note in passing that those operators per se are not 
strictly positive, but merely non-negative. Therefore suitable domain restrictions need to be carefully observed (e.g. 
this point is not a problem in the photon wave function construction, where p = is excluded by the transversality 
condition, see e.g. 0, H^l- 

Arguments of the present section are not necessarily limited by the dimensionality of space {D > 1 ) . Our discussion 
is based on the Fourier representation, hence we assume a priori the existence of (both direct and inverse) Fourier 
transforms of considered functions (the inverse transform is not a mathematically obvious notion, albeit always tacitly 
assumed to hold true in the phy sics literature) . We indicate how the domain of an inverse fractional Laplacian may 
be interpreted, see e.g. also |47| . 

The action of any | A|^/^ upon the Fourier transformable function f(x) can be defined as a multiplicative modification 
\k\^^'^f{k) of its Fourier transform f{k). Since the inverse Fourier transform [|A:|'^/^/(A:)]^ (x) = g{x) is assumed to 
exist, we have in hands an example of a "nice" function. For such function we can readily define the action of |A|~^/^ 
as \k\^'^^'^g{k) where g stands for the Fourier transform of a "nice" function g. In general, g cannot be arbitrarily 
chosen and rather severe restrictions need to be observed, to secure the existence of the inverse Fourier transform of 
|fc|^''/^g(A;). Our ultimate conclusion is that an identity |A|~^/^|A|'^/^ / = / surely makes sense for suitable functions 
f = f{^)- 

We are now ready to address the meaning of V ^|A|'^/^ necessarily appearing in the transformation of dtp = 
— lAI'^/^/o into a continuity equation dtp = — Vj. We consider k^{x,t) = p^{x,t), see Eq. ([57)) . For clarity of 
discussion we pass to the space dimension D = I. Accordingly 

■ r+00 

W^^lAl^"^^ Pf,{x,t) = — sgn{k)\k\i'-'^ exp{-ikx-\k\''X)dk (51) 

where the sign factor sgn{k) is a crucial ingredient that makes a difference between V~^|A|^/^ and | A|'^'^~^^/^. 

We can be sure of the existence of the above Fourier resolution for j{x,t) only in the stability parameter range 
p G [1, 2). For p S (0, 1), the function dtp^{x, t) cannot be represented as a gradient function. Therefore, for p £ (0, 1), 
the transport equation 9tp^(x,t) = — A^/'^p^j(a;, t) cannot be given a functional form of the continuity equation. 
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E. Pdf transport equations induced by confining Levy semigroups 

In case of "free" semigroups, there is no asymptotic invariant densities so that simiUar to the Brownian case, c.f. 
dtp = +Ap (we disregard so far irrelevant coefficients), we deal with a sweeping motion. 

The situation changes drastically, if we pass to the confining regime of a jump-type process, with an invariant 
density p*. It is rather obvious, [H, HI, Ull , that by choosing pl^^ and then making a formal replacement of (—A) by 
|^|A'/2 everywhere in Eq. ([8]), followed by an adjustment of noise intensity parameters 13 —> A we end up with: 

dtp = X[-pV'\A\^^'{p-'/' ■ ) + p-'/\\Ar^'pV')] p . (52) 

We note that, instead of a specific choice jAj''/^, we can equally well employ another admissible F(ji), like e.g. 
V— A + m? — TO. Our arguments extend to those cases as well. 

This time evolution of the pdf p{x, t) is formally induced by the (would-be) Lcvy-Schrodinger semigroup exp(— iif^). 

At this point we further extend an applicability of our arguments by considering a fractional (pseudo-differential) 
generalization of the "free" Hamiltonian (1) (set e.g. —DA — > A|A|^/^ and V/2mD — > V) to that with a confining 
potential: 

H^ = X\A\^/^+Vix). (53) 
A tacitly presumed ground state pl^^ corresponds to its bottom eigenvalue 0, i.e. H^pl^^ = holds true. Accordingly, 

at*==-iJ^^':=:-A|A|'"/2*-V(a;)^', (54) 
where ^(x, t) = p(x, t) p* (x). A compatibility condition, 

IAIm/2 

V = -A l^^l J^* , (55) 

P*' 



determines V and that in turn allows to expect an equivalence between Eq. ((52|) and Eq. ([54| . 

There are however some jeopardies in this formal procedure. An independent check is necessary of whether H > 
actually is a well defined self-adjoint operator with a dense domain of definition, which may not be the case. 

Things are much easier, if we first choose an appropriate iJ, together with its bottom eigenvalue corresponding 

1/2 

to the (ground) eigenstate . Then the semigroup is given a priori. Consequently, Eq. (j52p follows along with a 
proper large-time asymptotic behavior of p{x,t) — > pt{x). 

One can rewrite the transport equation (j52p in the canonical form of so-called master equation, appropriate for 
jump- type processes, [H, [s^ : 

dtp = [w{x\z)p{z) — w{z\x)p{x)]dz . (56) 



If we replace the symmetric jump rate, like e.g. (ID stable case is exemplified for clarity of arguments) 

w{x\y)=wix\y)^l/\x~y\^+'' (57) 

by a non-symmetric expression 

exp[^(x) - <i>{y)] 
^-^(^ly) |^„^|i+M 

then |A|^/2 |A|(^^^ and the corresponding transport equation reads: 

{1/X)dtp = -|A|^/V = -cxp(<f>) |Ar/2[exp(-$)p] 

-t-pexp(-$)|A|^/2gxp($). (59) 

1 /2 

We can always always select exp($(x)) = pj in the above, so arriving at Eq. ([52]). 

Remark: Independently of the method adopted, an important difference needs to be spelled out, if compared with 
the diffusion-type reasoning. Namely, the transport equation Eq. (|52p cannot be converted to the Fokker-Planck form 

dtp^-w(-^p)-X\Ar/'p, (60) 



TO/3 

traditionally attributed in the physics literature to Levy flights in external conservative force fields {U{x) stands for 
the Newtonian potential, [ioMS ]. (Note a replacement -DA -)■ X\A\f^/^ in dtp = -y{bp) + DA p.) This specific 
dynamical inequivalence issue has been pointed out in [l^ll^ll^, see e.g. also [l^. ISTj. 
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IV. NONLOCAL QUANTUM DYNAMICS. 
A. Duality transformation 

Fractional Hamiltonians = 7|A|''/^ with < /i < 2 and 7 > and likewise the quasi-relativistic one Hm = 
\J + rri^ — m with to > are self-adjoint operators in suitable L?{R) domains. They are also nonnegative operators, 
so that the respective fractional semigroups are holomorphic (also named analytic), and we can replace the time 
parameter f > by a complex one a = t-\-is, Rea > 0. Accordingly, a holomorphic extension of the Levy-Schrodinger 
semigroup is defined as follows, c.f. Eq. (pi]): 

[exp(-a7?)]/ = [(/cxp(-aF)]^ = / * fc. ■ (61) 

Here, the integral kernel reads = [exp(— ctF)]^ . Since H is selfadjoint, the limit i J, leaves us with the unitary 
group exp{—isH), acting in the same way: [exp{—isH)]f = [fexp{—isF)]^, except that now kis ■= "y|^[e2^p(— *s-F)]^ 
no longer is a probability measure (transition probability density). 

In view of the unitarity, the unit ball in is an invariant of the dynamics. Hence probability densities, in the 
standard quantum mechanical form p = ip* ip can be associated with solutions of the free fractional (or Salpeter) 
Schrodinger-typc equations: 

idsij{x, s) = 7I A|^/2^(a;, s) (62) 



idsip{x, s) = [V— A + TO^ — to]?/;(x, s) (63) 

with initial data ijjjx, 0). Attempts towards formulating the so-called fractional quantum mechanics can be found in 
Refs. BSEllii. 

We note the nonlocal action of motion generators is somewhat blurred in the Fourier representation. The pertinent 
spatial nonlocality becomes obvious if the canonical quantization is carried out on the level of the Lcvy-Khintchine 
formula, c.f. subsections HI. A and III.B. 

Guided by Q we identify the semigroup time label t > with the Levy-Schrodinger time label, e.g. set s — )■ t. All 
that amounts to the duality (Euclidean) mapping t ^ it oi Section II. D., c. f. (1) and p5|) . which we exemplify for 
the symmetric stable noise generators: 

dtp = -7|A|'^/2p — ^ idtt/j = 7|A|^/V (64) 

Stable stochastic processes and their quantum counterparts arc plagued by a common disease: it is extremely hard, 
if possible at all, to produce insightful analytic solutions. To get a flavor of intricacies and technical subtleties involved, 
whose neglect leads to erroneous formulas (and a danger of untrustworthy physical conclusions), we have been quite 
detailed in the analysis of Levy dynamical semigroups and their unitary (quantum) partners. Subsequently, while 
developing a general theory, we shall pay special attention to the quasirelativistic (Salpeter) equation and its m — !• 
Cauchy-Schrodinger limit. For clarity of arguments (and computational convenience) unwanted parameters (like 7) 
will be scaled away. 



B. Fourier representation advantages and drawbacks 



1. Levy-stable case 

Equations (|32p and ([55]) define the action of Levy stable and quasirelativistic generators upon functions in their 
domain. The involved integrals are interpreted in terms of their Cauchy principal values. One must as well keep in 
mind a crucial role of the ~ ~f{x) counter term in Eqs. (31) through (34). Its presence there is indispensable and 
appears to have beeen overlooked in quantum mechanically oriented papers, c.f. [30l. Is^. 

To explain that issue, we shall first discuss in some detail a fractional < /i < 2 Laplacian in space dimension 
n = 1. Its spatially nonlocal action upon functions in a suitable Hilbcrtian domain (domain issues we relegate to the 
last section of the paper) reads: 

i^Ar^^m = I Ar/v(.) . i vr/(.) = _ ni+,)sH.,/2) r m-m 

J \x-y\'-+t^ 
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Let us investigate the properties of — | A|''/^/(a;) by turning over to the Fourier image f{k) of f{x). Eq. (|65|) yields 
(a formal interchange of integrations is here-by executed): 

The integral over dy, presuming its very existence (which is not the case for = 1), can be calculated as follows 

We note the importance of the restriction fj, G (0,2) and obvious divergence problems to be taken care of: the T 
function is known to have simple poles at points 0, —1 and —2. Therefore, at /i = 0, 1, 2, the integral ([S7)) is divergent. 
However, irrespective of how close to 0, 1 or 2 the label /i > is, the integral (p7|) is well defined. 
It is interesting to observe that the divergence of the the Fourier integral, as /i approaches 0, 1 or 2, becomes 
compensated, if we substitute it back to Eq. (j66|) and next consider the limiting behavior of the result: 

2r(l + w)r(-u) sin cos r°° 

n\/2TT J -OC 

\k\^'}{k)e'^''dk. (68) 



Here we employ the identity (Euler's reflection formula), with an obvious reservation that it becomes invalid at sharp 
values /i = 0, 1, 2 of /i e [0, 2], while being operational for all /i G (0, 1) 2): 

r(i + Ai)r(-M) = — (69) 

Sm TT/i 

This observation clearly identifies some of frequent misuses of the formalism if sufficient attention is not paid to 
potential obstacles (like e.g. an interchange of improper integrals or a neglect of the F function simple poles along the 
negative semi-axis). For example, the range of validity of the right-hand-side of Eq. goes beyond ji £ (0, 2) and 
admits a safe extension to the boundary values and 2, by-passing as well the previously raised problem concerning 
A^ = l- 

Quite apart from this appealing outcome, the primary integral representation ([65]) is not valid at the boundaries 
of the stability interval and one needs to resort to the remaining local terms of the general Levy-Khintchine formula. 
pT]). More than that, if one presumes the Fourier representation Eq. (|134p as a valid definition of how |A|'^/^ acts 
upon functions in its domain, there is now clean way to go backwards, such that the primordial definition (|65l) actually 
could have been reproduced. 

This obstacle has been often overlooked in the literature. That can be explicitly seen in publications on the 
Salpeter equation and its solutions, [13, HI]- For example, in the would be (actually divergent) integral kernel 
formulas for operators |A|'*/^ and [V— A -I- — m] respectively, presented there-in, the ~ ~f{x) counter term is 
conspicuously missing. Below we shall be more explicit on this " missing counter term" issue in the discussion of the 
ID quasirelativistic case with to > 0. 

2. Quasirelativistic Schrodmger (Salpeter) equation and its m ^ limit. 

In the physical units, the ID relativistic (here named quasirelativistic) free Schrodinger equation is commonly 
considered in the form 



Denoting /(fc) = (27r)~-^/^ /(a;)e~''^'^dx, f{x) — (27r)^^/^ /(fc)e*'^'^dfc and interpreting the action of the 



square root operator ymPc^ — ft^c^-^^ '/'(s;,^) in terms of the series expansion, we readily arrive at the following 
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formal Fourier representation: 



92 



/27r J- 



(i>{k,t)dk 



1 92 



mc 



V2^ J- 



<j,{k,t)dk Wl + ^ e'''- 



m2c2 2 9a;2 
1 



1 92 



ikx 



ikx 



(71) 



We note that h/mc stands for so-called reduced Compton wave-length and the momentum label has physical dimen- 
sions p = hk. 

Although we have anticipated the existence of the mass m = limit in the relativistic Hamiltonian context, the 
above derivation of (j7ip rings warning bells. Indeed, tacitly presuming the nonrelativistic regime m2c2 we have 
plainly expanded toc2-^1 + [p'^ /nrfic?) in ([71]) into Taylor series with respect to ~ p^ /m?c^ and evidently we are left 
with no room for m — > therein. 

Nonetheless, we can safely put m = 0, after the series resummation - in the last entry of the formula (j7ip - so arriving 
at the correct form of the Fourier image of |V|. Indeed, [l^l, to this end we should consider the ultrarelativistic regime 
with p2 ^ 77j2 and make an expansion of the \p\c ^Jl + (w? (? j p^') with respect to ~ rr?(? jp^ . Letting m ^ becomes 
a legitimate operation that replaces \/ m^d^ ~V p^c^ by c\p\ in the Fourier representation. We refer to the previous 
subsection for a discussion of how c|V| can in turn be recovered, c.f. also Eq. (78) below. 

We point out that the nonrelatvistic limit does make sense exclusively in the Fourier representation. More than that, 
the limit m — >■ and Fourier imaging of Eq. (|70p are not interchangeable operations. This is an important subtlety of 
the mass to > quasirelativistic dynamics, quite akin to those raised in relation to Levy stable Hamiltonians of the 
previous subsection. 



C. On integral (kernel) representations of quasirelativistic and Cauchy generators. 

In the present subsection we take under scrutiny a procedure [sol . [3^ of assigning an integral kernel [not a 



semigroup kernel discussed in Section III) to the operator im? — and likewise to fi,c|V|. An origi- 

nal argument of Ref. [lOl goes as follows: (i) assume ([71]) to hold true, (ii) take the inverse Fourier transform 
^(p,t) = (27r)~-^/2J' (j)(^y^t) ex]){—ipx/h) dy, (iii) use an identity \J + x"^ exp{±ipx)dx = —{2a/\p\)Ki{a\p\), 

The outcome presented in Eq. (2.24) of Ref. HO], see also Eq. (14) in Ref. [11], reads ^Jra^c"^ — h^c^-g^ 4>{x-,t) = 

K{x — y)(f>{y,t) dy. This result is plainly incompatible with the primordial formulas p3p and (|43p . The same 
comment (concerning the faulty outcome) refers to an analogous reasoning for /ic|V|, with the ultimate kernel K{x — 
y) = —chj'Kix — j/)2. The crux is that the above arguments have not been properly worked out. 

Namely, while departing from the above mentioned identity (iii), let us evaluate an auxiliary integral: 

^^[°|^'^ (e'"P - \)dp = — r dp r dx e^P-^/^^^^ie'^P - 1) = 
\P\ 2™ J_^ 

a2 dx / (e*P(^+'=) - e'f^) dp = - [S{x + k) - 5{x)] yGFT^dx = 




2 



a 



(72) 



Since we actually have 



92 TTjr2 ^ / r)2 

n2c^ ^ ^(a;, t) = ^ 0(fc, t) . / 1 + ^ e^k. (73^ 
ox^ V27r J_oo V ™ c2 
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and an identification a = mc/h implies ik/a) = p /m c\ we realize that 



7rV2^J-oo \V\ 



= mc'cf>{x,t)^— I + (74) 

Ti" i-oo |y| 

This implies a correct integral form of the ID Salpeter equation for 'ip{x,t) = exp(imc^t//i) <): 



zft — = U / m^c* - /i^c^ — ^ - TOc ?/>(a;, i) = / [ijj[x + y,t)- -iplx, t)\dy (75) 

dt V dx^ TT |y| 

Compare e.g. Eqs. (gS]), see also Eq. (15) in Ref. 

Some comments are in order here. First, we pay attention that in fact we have accounted for the "missing coun- 
terterm" from the very beginning of our calculations, by inserting a factor (e**^^ — 1) in the integrand. One should 
realize that, if taken literally, the integral 

J-oo \y\ 

is divergent, even in the sense of the Cauchy principal value. This is a consequence of its l/p^ behavior of the integrand 
as p — >■ 0, and its evenness. On the other hand (set a = mc/h) the corresponding "regularized" integrals ([72]) and 
([7i)) are finite 



1 rMiW ,,^^^74: 

TT \y\ V m^c- 



(76) 



For reference purposes, we list one more useful integral for ^ £ (0, 2), c.f. Eq. f|67p : 



1 poo 1 

\kr = 7^^^, WW / - !)• (77) 

2r(-/z)cos^ 7„oo 1^1+^ 



The Cauchy case is introduced through the limiting procedure ^ ^ 1, where r(— ^) cos —tt/2 and it follows 



(78) 



At this point we realize that the to — > limit can be executed in (74). To this end we must multiply (|76p by mc^ 
and take notice of Ki[x 0) — 1/x. Accounting for p = hk ultimately reproduces ([75)). next see Eq. (|134p . The 
corresponding Cauchy-Schrodinger equation reads 

J-oo \y\ 

and provides a consistent definition of a spatially nonlocal generator of quantum dynamics (Eq. (|127p likewise). Our 
discussion, including that of section IV. B, clearly demonstrates that the Fourier transcription of spatially nonlocal 
expressions (like e.g. Eq. is a derived secondary ingredient of the theory which must not be employed hastily, 

but with due care. 



D. Propagators in ID. 



Levy semigroup kernels have the general form Eq. ([57)1 and ([551 The duality transformation t ^ it replaces the 
semigroup operator by an affiliated unitary operator and likewise the Fourier representation of the the semigroup 
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kernel by that of an affiliated propagator (transition amplitude, unitary group kernel). Accordingly, keeping in mind 
that in general k = k(x — y,t — s) with t > s > 0, we have: 



k^ix,t) kf,{x,tt) = K^{x,t) = (27r)-" / exp(-ipa; - ii|p|^) d> 



(80) 



and 



km{x,t) km{x,it) = Km{x,t) = (27r)" 



exp[ 



-ipx — it{\/ m? + — m)] (Pp . 



(81) 



see e.g. [131 fo'^ ^ parallel discussion of the Gaussian (heat kernel vs free propagator) case. 

First of all we must take under scrutiny, quite often met in the literature, appealing but naive short-cut that 
amounts to mimicking the t ^ it substitution of Section II. D directly in the spatial expression for the Levy semigroup 
kernel. This procedure could have been justified in the Gaussian case [l3| but is invalid if extended to Levy kernels 
without suitable precautions. 

The problem is that a formal analytic continuation in time of inverse polynomial pdfs typically produces singular 
functions with poles. Time honored quantum field theory procedures were developed in the past to handle similar 
pole problems, but no mention nor trace of them could have been found in a number of papers devoted to relativistic 
quantum mechanics, [21_, 28, 3Q, [s^ . 

A fairly typical example of repeatedly reproduced erroneous outcomes are e.g. formulas (4.31) and (4.34) of Ref. 
[2ll | and likewise (29), (C.5) and (C.19) in [l^. In particular, for the would-be ID Cauchy propagator, the faulty 
formula t/Tr{x'^ + t^) —J- (l/27r) J cxp(±ipx — it\p\) dp = it/TT{x'^ — t^) is reproduced. We have addressed the involved 



Fourier integral before, [8[. Presently, we shall give two complementary derivations, independent from the previous 
one, both resolving the poles problem. 

We take the previously mentioned Cauchy kernel kt = tj-Kix^ + 1^) and rewrite it in the form 



ktix) 



I 

2^ 



1 



1 



X + it X — it 



(82) 



To arrive at the Cauchy propagator, we perform a formal substitution i it + e, to be followed by the limit e — > 0. 
This yields 



Kt{e) 



I 

2^ 



1 



1 



X — t + is X + t — ie 



(83) 



In view of the well known identity, |45j . (P{l/x) indicates that the generalized function 1/x needs to be interpreted 
in terms of the Cauchy principal value of the involved integral) , 



1 



xztie 



= P 



^ i'!TS{x), 



(84) 



we get 



Kt{x) = - 
Zn 



P 



1 



^P 
2n 



X — t 
1 



- iTTS{x -t)- P 



1 



X + s 



x + t 
6{x — s) + S{x + s) 



— i'K5{x + s) 



= -P 

TT 



t2 



S{x-t) + S{x + t) 



(85) 



which is a correct expression for the Cauchy pr opagat or, previously obtained in Q. That needs to be compared with 
the naive outcome it/TT{x'^ reproduced in [2l|, l28|. We note that the e-regularization still survives in the formula 

(|85p . although is not explicit. Its tacit presence is somewhat blurred by the Cauchy principal value indication and 
the emergent Dirac deltas. 

Let us represent the e-rcgularized Cauchy kernel Eq. (|83|) as the following Fourier integral 



1 

2^ 
1 

2^ 



OO 
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2^ 



OO 1 

e~'^'~"'^^''^ cosux du= — 

-OO 



27r 



-i(t-ie)\u\^±iux 



i{t — ie) 



n[x'^ — [t — ie)^ 



(86) 
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It can be seen from ([86]) that the e - regularization is introduced to secure he convergence of integrals quantifying 
the wave packet evolution by means of kernel functions. After an explicit evaluation of the integrals, we can safely 
put e — > 0. In view of this implicit e — > limit, the — ie term in the numerator is in fact irrelevant and can be safely 
neglected, yielding a familiar form of the regularized ID quantum Cauchy propagator 

^* = ^[x^~c^{t-ieY] ■ ^^^^ 

where we have reintroduced the velocity of light c-dependence. 

In the literature, one encounters an appealing (in view of the discussion of Section II. D) but faulty mapping of 
the Cauchy transition density into the Cauchy propagator kt = t/'K{x^ + 1^) Kt = it/n{x'^ — t^), with no mention 
of the pole problem and the need for a regularization, see for example [U, In fact, if we tacitly disregard the 
pole obstacle and, while on the level of the Fourier representation, formally set e ^ in K^{x) and next perform 
integrations in Eq. ([50)) (specialized to n = /i = 1), the outcome would actually be Kt = it/Tr{x'^ — t^). As indicated 
above, the fully-hedged duality transformation exp(— <|V|) — >■ exp(— it|V|) enforces an e-regularization of the kernel 
function, implying less straightforward, but undoubtedly correct outcome Eq. (|86p . 

The above reasoning extends to the quasirelativistic case as well. Namely, for the pertinent semigroup 
exp[— t(V— c^A + ni^d^ — mc^)] the ID kernel function (transition probability density of the associated jump-type 
process) reads: 

fc™(x, t) = cxp(TOc2i) ^ (88) 

TT ^Jx^ + cH^ 

By turning over to the duality transformation t — > it, while remembering that K\(z) ^ 1/ z for small values of 2, we 
can readily produce a regularized version of the quantum propagator (note that we use = again) 



^ , , m(?t , T , K\{mc\/ x^ — c?{t — ie)^) 
K'(x, t) EE i exp mc^it ^ , ^ ^ ' ' . 89 

whose m — > limit coincides, as should be the case, with the previously defined quantum Cauchy propagator. By 
turning back to the semigroup transport formulas (j34p , (j35p and I36p and directly executing the duality mapping 
t ^ it right there, we can easily discover an origin of the e-regularization, that is necessary if one attempts to deduce 
a consistent analytic expression for quasirelativistic and Cauchy quantum propagators. Let us focus on the Cauchy 
case again. The unitary dynamics of any ip^x, 0) = f{x) is obtained as follows 

r+oo 

ip{x,t) = [cxp{-it\V\)f]{x) = (271)-^/^ / exp{-it\k\+ikx)f{k)dk = {f'>'Kt){x) (90) 



where (/ * Kt){x) :— J Kt{x — z)f{z)dz and Kt — (l/27r) J cxp{—ikx — it\k\)dk. As long as we do not insist 
on an explicit evaluation of the Fourier integral expression for Kt{x), no regularization is necessary. The Fourier 
expression for 'ip{x,t) evidently does the job without any special precautions. On the other hand, would we have 
turned to ^{x,t) = (/ * Kt){x) := J Kt{x — z)f{z)dz, it is only the regularized expression (jSS]) that yields a time 
development of il>{x,t), in conformity with the Cauchy-Schrodingcr equation, The naive would-be propagator 
^najtie^^j _ j^^^^^2 — t^'^^ if taken literally, docs not lead to a consistent evolution pattern, see below. 



E. Quantum wave packet dynamics in ID 

1. Heavy-tailed initial data 

With the propagator in hands, we can address the dynamical behavior of solutions of the free fractional (eventually 
Salpeter) Schrodinger-type equation. Lets us continue a discussion of the specific ID Cauchy-Schrodinger case 

idti!{x, t) ^ I A|i/2V'(a;, t) = |V|V'(a;, t) (91) 

As an initial condition we take an "almost Lorentz" (up to normalization) distribution, so that |?/'|2(a;,0) actually is 
a normalized quadratic Cauchy pdf, [39| : 

V'(x,0) = /(x) = ./^-^. (92) 

V TT 1 + 
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FIG. 1: Time evolution oi p{x,t) — \ip{x,t)f, (I95|l . for the Cauchy-Schrodinger dynamics (|9ip with the initial data (|92p . Figures 
near curves show time t values. 



The integral kernel of exp(— i|V|t) takes '^(a:, 0) into ip{x,t) through a convolution, c.f. Eq. (|36p . For clarity of 
discussion we employ a regularized form K^(x) of the integral kernel, allowing to control its behavior in the vicinity 
of singularities, so arriving at 

fix - z)Kt[z)dz =Tr\- T—i ^ / cosuz du. 



Consider I ^ JT^ ffp^ = ^^^^^^^ 



and an explicit evaluation of the integral, followed by e — > 0, gives rise to 

^{x,t) = + t) + fix - t)] t)f{x ^t)-{x + t)f{x + t)] = y|^^^d_iL_, (94) 

where f{x) stands for the initial data for the pertinent evolution. Taking the modulus of a complex function ([M)) 

1/2 

we get interesting outcome, with an implicit redefinition {x) = f{x): 



\^{x,t)\- = (1 + t-Wp,{x + s)p,[x - - - ^ , ^ ^ , (95) 

which conveys a useful information about an enhanced spreading of the initial wave packet due to the dynamically 
generated bimodality of the square root expression. Compare e.g. [1, [2^ HO, HH . The pertinent time-evolution is 
visualized in Fig.l, compare e.g. also [28| . 

We would like to point out a definite advantage of the (e) regularized Fourier representation of the kernel which 
permits to incorporate naturally the delta-function terms and simplifies calculations. A complementary detailed 
derivation, explicitly accounting for delta-type contributions (j85p . can be found in Ref. Q. 
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2. Link with a classical (d'Alemhert) wave equation. 



The solution (|94)) of Eq. ((9T|) at the first glance seems to look like a general solution of the original equation 
of motion, represented as a superposition of special solutions (wave packets that are "moving" right or left in the 
standard quantum lore), like e.g. <^{x,t) ~ {l/\/2[(j){x — t) + 'iIj{x + t)], c.f. Eq. (3.37) in [s^l- This is merely an 
illusion and wc must consider (|94p as an indivisible whole, since no physical status can be assigned to separate right 
or left "moving" components. This is a direct consequence of the manifest nonlocality of the motion generator |V|. 

Our standpoint is that there is no such entity as the right or left "moving particle" in the present framework. That 
stays in plain opposition with a discussion of Sections 3 and 4 in Ref. [30|- Even worse, the very "particle" notion 
(concept) appears to be doubtful in this nonlocal dynamics setting, in view of the well known "particle" (especially 
massless "particle") localization problems plaguing the relativistic quantum mechanics, c.f. for example Section V of 
Ref. [Ijand 0,1131 and references there- in. 

Eqs (j94p and (j95p refer to a nonlocal complex-valued wave phenomenon, implying that the induced probability 
distribution bimodally expands and ultimately fades away (spreads). The underlying dynamical mechanism amounts 
to the propagation of local maxima, in the opposite directions, with the velocity of light (here c = 1), ([55]) . None of 
those maxima separately can be given a status of physical relevance, e.g. must not be interpreted as an identifier of a 
right or left "moving particle" wave packet. 

On the other hand the left- and right "moving" components of (|M)) can be given a physical status within the higher 
level theory, i. e. the pure wave ID d'Alcmbcrt equation (— A-(- ^)u(x, t) = 0, where c = 1. Both pertinent "moving" 
components are solutions of this wave equation and likewise their superposition (j94p is. The crux is that neither of 
component functions is by itself a solution of the Cauchy-Schrodinger equation. It is the "indivisible whole", i.e. 
their superposition ipi^jt) of Eq. (|M)) . which for sure solves ([?T|) and admits a standard quantum mechanical (Born's) 
probabilistic interpretation. 

It is worth pointing out that the Hilbertian domain of the operator |V| (containing all normalizable solutions of 
dtt/j = I Vji/') can be consistently built by putting through a suitable sieve the set of all solutions of the ID D'Alembert 
equation. The pertinent domain is a fairly restrictively selected subset (closed linear space) of solutions of the 
d'Alemebcrt equation, such that the Cauchy-Schrodinger equation is solved by them as well (that we know not 
necessarily to be the case). 

Remark: With Eq. (|94p in hands we can be more explicit on the last sentence of the previous subsection. 
Namely, the Cauchy principal value of the (otherwise divergent) integral representing a "naive" propagation ^(x,t) = 
(/*i^""™*^)(a;) := / K^^™'^{x — z)f{z)dz of 4'{x, 0) = f{x) = ^J2/^T / {l + x"^) would result merely in the pure imaginary 
term [i / 2)[{x ~ t) f {x — t) — {x + t)f{x + t)] of the expression (|94p . which we know not to solve the Cauchy-Schrodinger 
equation as it stands, c.f. |8| for more details. 

3. Dynamically generated bimodality 

In view of Eqs. (|94l) and ((95|) . a minor modification of the initial data ([92|) to ^p{x.O) = •^^^[7/7^ -f- x^], implies 
^p{x, t) = [(7 + *0/[(7 + **)^ + ^^]] and thence gives rise to the time-dependent pdf of the form 

" + 2^2(72^- i2) + (^2+^2)2 ■ (96) 

We can quantify an emergence of bimodality by investigating an extremum of p{x, t) with respect to x, i. e. solving 
an equation dxp = 0. This amounts to solving x^ + ^(7^ ^ t^) ^ with an obvious outcome: xq = 0, and two more 
real roots defined by x± = ±y^t'^ ^7^ under the condition t"^ — > 0. 

Accordingly, for times t^ < 7^ the considered pdf is unimodal, while at time instants — _ q ^-^^ situation 
changes. The bimodal form od the pdf is born and persists for all t^ > 7^. The pdf of the previous subsection 
is recovered by putting 7 = 1. Wc shall show subsequently that the reported bchaviot of Cauhcy-Schrodinger 
pdfs is not due to a special choice of initial data. Additionally, we point out that the dynamical generation of bi- 
modality is not special to the Cauchy-Schrodinger case and appears as well in the quasirclativistic (Salpeter) evolution. 

Remark 1: There is no qualitative change in the behavior of solutions if we pass to the ID quasirclativistic 
(Salpeter) equation dttp{x,t) = V— A + ni?'ip{x, t). Initial data of the form 



^(x,0) 



7rA'i(2m7) 



7^:1 (m^a:^ -f 7^) 

AT"; — 2 ^ ' 
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FIG. 2: Time evolution of p(x,t) = \ijj{x,t)\'^ for the Cauchy-Schrodinger dynamics (|100|) . with the initial data (|98|) . Figures 
near curves show time t values. 

evolve in time according to a simple substitution rule 7 — >■ 7 + zt. As well, one easily verifies that the m — > 
limit reproduces the Cauchy-Schrodinger wave function, whose is indeed the pdf p{x,t) of Eq. (|96)) . see e.g. 

[28I [30I Issj . The denominator -^/a;^ + (7 + it)^ is a source of an emergent bimodality of this particular solution, as 
first noticed in Ref [33j . 

Remark 2: For both ID Salpetcr (m > and m ~ 0) cases, a continuity equation dtp = — Vj has been verified to 
make sense and suitable probability currents arc known in a closed analytic form, see e.g. [soj . 

4- Gaussian initial data 

Now we consider the Cauchy-Schrodinger time evolution of the initial Gaussian wave packet, c.f. [2^ for a comple- 
mentary discussion of this issue, 

^(x,0) = ^e--V4. (98) 
The dynamics is made explicit by employing the Fourier representation of ip(x^t). Namely, we have . 

V2^y-oo V2^V^7-oo 

9 \ 3/4 „oo 



cos(za;) dz. (99) 



The integral can be evaluated exactly in terms of special functions to yield 



where crfi(x) — —i erf(ix) is a real function. 



e'^^crfi ( I + crfi ^' ^ ^ 



(100) 
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FIG. 3: The probability current for an initially Gaussian distribution in the ID Cauchy-Schrodinger dynamics. Figures near 
curves show the time t values. 

We note in passing that following analytic methods of Ref. Q we would obtain another explicit form of the above 
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(101) 



The squared modulus (e.g. a corresponding pdf) of the expression (|100l) (or (|10ip ) is displayed in in Fig. 2. A 
qualitative similarity between Gaussian and Cauchy (|95p cases is transparent. Both functions exhibit a dynamically 
generated bimodality, see also [l^. This behavior oi ip{x, t) is definitely dictated by the propagator and appears not 
to depend on the initial data ip{x^ 0) choice. 

We finally note that there is no problem to deduce, with a numerical assistance if necessary, the detailed time 
evolution of any initial wave packet (Gaussian, Lorentzian etc) driven by a quasirelativistic kernel (e.g. proportional 
to the MacDonald function Ki). 

Let us calculate the probability current for the present case (e.g. Gaussian initial data for the Cauchy-Schrodinger 
dynamics. To this end we shall employ the observations of [s^l. The Fourier transform of ^'(a;, 0) = "i^^^^^ reads 

'0(^,0) = ^J^e~^ . The evolution can be encoded on the Fourier transform level as follows: '0(p, t) = ^J^e~^ -*b|t. 

Next, we have to insert ip*{p,t) and ipik^t) to the defining identity (see e.g. subsection IV. E, here we scale away all 
dimensional constants and consider ID instead of 3D): 



1 

2^ 



dp 



dkJ^e^^ 



'\p\ + \k\ 

To facilitate the calculation of integral we pass to the polar coordinates p = rcosLp, k — r sin 1^9. Accordingly 

1 



7rv27r 



Sm ip + cos If f°° ^^-r^ ^ixr{sinv-cosv)^itr{\cosip\-\sin<p\)^^^ 

I sin (/3 1 -I- I COS V3 1 Jq 



(102) 



(103) 
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The integral in (|103p can be evaluated analytically 

j- , X{(p) = X {sin (fi — cos (p) + t{\ cos (p\ ~ \smip\) , (104) 

so that the final expression for the probability current assumes the form (cross-checked by means of Mathematica 
routines): 

^ 1 f^^ smtp + cosip 

^^'''^"(2^X |sin^| + |cos^| 

A numerical visualization of the time development of that current is reported in Fig. 3. It is seen that behavior of 
the current is qualitatively similar to that for the Cauchy initial packet, see e.g. Fig. 2 in Ref. [s^l- Namely, at small 
times, there is practically no current at all (we begin from j{x,t = 0) = 0). The current is an odd function of x. As 
time passes, two positive and negative current peaks develop and move towards plus and minus infinities respectively, 
while falling down to zero. 



re 



dr 



crfi 



1 e 



dLp. 



(105) 



F. Quantum wave packet dynamics in 3D. 



1. Generic wave packet evolution pattern: radial expansion. 



Peculiarities of the 3D dynamical behavior can be conspicuously seen in the Cauchy dynamics idt'tp = |V|'!/', where 
due to the radial symmetry, the ID bimodality is replaced by a "spherical" -modality (and "circular" -modality in 2D). 
The initial wave packet rapidly delocalizes by a dynamically developed expansion of a spherically-shaped probability 
concentration to the radial infinity, mimicking an expanding spherical wave with the source at the origin. Its radial 
peak is running away with the velocity c (here interpreted as c = 1): the pdf radial maximum resides on a surface of 
an expanding sphere (an expanding circle in 2D and an expanding interval in ID). 

This property can be directly deduced from the analytic formula (?■ = -^/x^ + + z^) 

implying (set 7=1) the following functional form of the related pdf: 

p(r,0) = |^(r,0)P = ^(YT^ = l^^^'^)"' = ^IP^TO^TW ^'"'^ 

Its time evolution is depicted in Fig. 4. (7=1,7" = ^y x"^ + y"^ + Z'^) . 

We note that the normalization coefficient reflects the radial symmetry of the problem (47r factor is the remnant of 
angular integrals), and Fig. 4 actually depicts the behavior of r^p(r, t), instead of the "plain" p(r, t). That should be 
compared with the time evolution pattern reported in Fig. 9 of [3 Of . 

The probability current for this case has received a closed analytic form in (30j . Both the dynamical behavior and 
the rapid decay of j{x,t) are depicted in Fig. 10 of [soj . 

Remark: Since a solution of the Cauchy-Schrodinger equation is also a solution of the D'Alembert equation (and 
not necessarily in reverse, c.f. subsection IV. E. 2), it is worth pointing out that the pertinent wave equation is known 
to admit radially expanding (or contracting) spherical waves that originate from (or sink in) a point source (or sink). 
Superpositions of expanding and contracting waves are legitimate solutions as well. 



2. Propagators 



In analogy to considerations of subsections III.C and IV. D, we may implement the quantum propagation ip{r, 0) — > 
'0(r, t) by means of the 3D propagator. To this end let us recall that the Cauchy kernel (transition probability density 
of the Cauchy jump-type process) stems form the general formula (j4ip . specialized to n = 3. Accordingly, while 
making explicit the c-dependence, we depart from 

fco(x - X', t - t') = ^2[^^_^'~^2\t^t,)2^2 = / d'kcM~M^ - X') - \k\c{t - t')] (108) 
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FIG. 4: 3D Caudiy-Schrodinger evolution: spherical expansion, exemplified for the radial pdf r^p{r,t) 



where |k| = <^k^ + ky + k"^. For clarity of further discussion we shall keep in mind a space-time homogeneity of the 

corresponding random motion and retain x and t labels, instead of x — x') and t — t' respectively. 

The duality transformation exp(— c|V|t) cxp(— ic|V|i) can be accomplished on the level of integral kernels, 
provided we take care of potentially dangerous singularities (poles). Like before that can be secured by considering 
an e-regularized Cauchy-Schrodinger propagator in the form (compare e.g. also (s?!): 



(2^)3 



(Pk exp[— ikx — ic|k| (< — ie)] 



ict 



c2(i-ze)]2 



(109) 



where r = |x|. C. f. also our comment following Eq. (|86|) in regard to = versus — in the above formula. 

Perhaps it is instructive to mention (and make a direct comparison with Eq. (jl09p the so called complexified 
fundamental solution of the d' Alembert equation, [23| , which has the form (we set e instead of the originally employed 
real constant a) 



1 



(110) 



For e 7^ this function is analytic in the whole space-time. Would we proceed formally (it is not quite legal step in 
view of our discussion of the poles issue), and set e = 0, one arrives at the Hadamard fundamental solution which is 
singular on the light cone. 

Armed with the lesson of subsection III.D and the above (|108l) to (|109p realization of the duality mapping, we know 
how to handle the arising singularity problems in the 3D (m > 0) Salpeter case as well. Namely, the quasirelativistic 
jump-type process transition kernel (e.g. probability density, whose space-time homogeneity is implicit), in view of 
Eq. (gni) reads 



km (x, t) 
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— -r- / d'^fc exp[— ikx — t{-\/ c^k^ + m?c 
ttY J 



27r2 



exp(TOc^t) 



K2 {mcy/r'^ + cH) 



J.2 _^ ^2^2 



(111) 



where K2{z) \/z^ for small values of z. 
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Accordingly, the e- regularized propagator of the 3D Salpeter equation dtip = (\/— c^A + mc? — mc^)?/' reads: 

K2(rac\/ — (?(t ~ it)) 



rr?e' 
27r2 



exp(mc it) 



(112) 



where r ~ \/ x"^ + y'^ + z'^ and h = 1 has been presumed. One may readily verify the validity of the m — > limit, Eq, 
(ITO . 



5. Probability currents m 3D 



For the completeness of exposition we find instructive to present a detailed derivation of the probability current for 
the 3D Salpeter (m > 0) equation, whose main idea we borrow from Ref. [soj . The transition to a simpler ID case is 
obvious. We represent the Salpeter dynamics with an additive perturbation by an external potential 



Accounting for 



(113) 



(114) 



and its complex conjugate we can evaluate the time derivative of the probability density p ~ according to 
mO^^Ml+^^Msothat 



5^ 
dt 



(115) 



Our aim is to transform (jllSp into a regular 3D continuity equation — §f= Vj, p = Taking the Fourier 

transforms of each if)- entry in the right-hand-side of (jllSp separately, and presuming that the r-h-s actually stands 
for Vj, we get 



Vj 



ic 



(27r)3ft7 



d^p d^k(t)l(t)k 



y/rn^c^ + fc2 -|- y/ m?c^ + 



(116) 
(117) 



In view of (k + p)(k — p) = k"^ — p^ , we readily retrieve the functional (spectral) expression for j that obeys the 
continuity equation: 



d^p d^k- 



^f(k-P) 



{2TTfhf> J ^ Vm2c2 + k'^ + ^m-^c? + p^ 
The massless case arises if we set m = in pisp . The probability current for the massless case reads 

k-hp 



(27r)3;i6 



d'-p d-'k- 



Oke 



if(k-p) 



(118) 



(119) 



To have a direct comparison with some of our previous discussions, one should set h = 1 = c whenever necessary. 

Remark: In the present subsection we have discussed the quasirelativistic generator that is additively perturbed 
by a function (like e.g. a Newtonian potential) that secures a bounded-ness of the Hamiltonian in (|113l) from below 
and (via an additive renormalization) moves its spectrum to with as its lowest eigenvalue. A number of typical 
(e.g . simplest) spectral problems have been addressed in the past both for the 3D and ID quasirelativistic cases, 
[3ll I58l - l6l| and references there-in. Some special spectral problems (harmonic potential and the infinite well case) for 
the m = 0, were considered e.g. in [3l|, [H, [H, 13 [H, [g^ and f63l - [65j . There are open problems even in those simple 
cases; controversies concerning an analytic form of eigenfunctions and the ground state function in particular have 
not been resolved for a " particle in a box" , [4^ [gJ, |65| . Mathematical arguments [6^ , indicate that none of proposed 
so far solutions is correct, as confirmed in the very recent Ref. (s^ . 
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4- Non-existence of the continuity equation for the fractional quantum dynamics. 

Ill the standard quantum mechanics, where the minus Laplacian is recognized to serve also as the Wiener noise 
generator, the time derivative of p is routinely represented in terms of the continuity equation. The argument goes 
as follows: consider a normalized solution of the Schrodinger equation (we set here-by h/2m ~ 1), then 9t|V'P = 

z[V'*(v2)V' - vv^vi = - ipvv) = -vj. 

If we consider the /i-stable generator instead of the minus Laplacian, a consequence of the Levy-Schrodinger equation 
is: dt\4'\^ = — i(?/'*|A|''/^V — V'l^l^'^^V'*)- This equation readily applies to the |V| = |A|^/^ and extends to the Salpetcr 
generator A + m^, whose mass m = version the Cauchy generator actually is. The just derived 3D version of 
the probability current can be associated with those two special cases. 

The situation appears more troublesome when general Levy-stable generators |A|^/-^ are involved. We note that in 
this case there holds (we refer to 3D, dimensional constants are scaled away) 

_ ah/f _ l_ L*IA|M/2,/,_,„IAI''/2,/,*l _ ' [ ^3„^3, |fc|'^- bl'^ W,*^..'f(k-P 



dt h 



(2^)3 J ' \kY + W 



and there is no way to extract the crucial factor k — p, needed to identify (jl20p as a legitimate (spectral) expression 
for Vj , with a notable exception of /i = 1. 

A candidate expression for the probability current, albeit restricted to the stability parameter range /i £ (1,2), has 
been produced in Ref. ([H, E^, see also an explicitly uncompleted passage from (38) to (39) in [H^- While setting 
h=l and disregarding the noise intensity parameter I?^, j has been presented in the form: 

j = -i?A*|A|''/2-iv?/. - VIAP'/^-iy?/.* (121) 

The problem is that the identity dtp = — Vj is invalid as it stands, in its differential form, e.g. there is no —divj local 
representation of dtp. This can be demonstrated by inspection after passing to the spectral representation. 
Our negative argument goes as follows, [H^. First we observe that by taking the divergence of j, (|12ip . we get 

Vj = -i [(VV-*) |A|^/2-iy^, _ (y^) lAI^/^-^W*] - i VIAI^'/^^VV- - i'\/\A\^'/'^-^\/^p*] (122) 

Now, we shall consider the second square-bracketed term in the above. To this end let us explicitly evaluate the action 
of involved operators, while passing to the spectral (Fourier) representation of the relevant entry. We have e.g. 

^*v|Ar/2-ivv = -^*(x)^^^ J d^p\pre^p-i^ip) = -r\Ar^^^ im 

and taking into account a complex conjugate of (|123p . we recover the right- hand-side of the transport equation (|120p . 
(up to h = 1 and an obvious sign inversion). For j of Eq. (jl2ip to be a legitimate probability current we must 
guarantee that the first square bracketed term of Eq. (|122p identically vanishes. This is not the case for complex 
functions ip. 

Remark: It has been demonstrated in @, c.f. also @, that the transport equation = 
-jr [7/.*|A|^/2^-l/'|A|''/2^*] can be literally represented as a master-type equation for the jump-type process with 
a corresponding jumping rate q(pi,t,A) where A is a Borel set in R^. The pertinent master- type equation has the 
form dtp{A,t) = J q{y,t, A) p{y,t) d^y and q{x,t,A) critically depends on the involved Levy measure. If p{x,t) is 
a probability density of a certain stochastic process, then p{x,t)d'^x = p{A,t) tells us what is the probability for 
a jump to have its spatial direction and size matching x G A. However, one should not be inclined to think that 
a concrete jump refers to a "physical particle" that jumps in space. We shall leave this debatable point open for 
further discussion. To the contrary, in a conventional Laplacian-based quantum mechanics, a standard interpretation 
of p{x,t){Ax)^ as a probability to locate a "physical particle" in a cube of volume (Ax)^ seems to be consistent. 

V. DOMAINS FOR NONLOCAL OPERATORS, STATE VECTORS AND UNCERTAINTY RELATIONS 

A. Mathematical prerequisites 



Hamiltonian operators discussed in this paper are unbounded. When defining an unbounded operator, it always is 
necessary to specify its domain of definition. We proceed in he following general manner. If A is an operator in the 
Hilbert space H, we write D{A) C H for the domain of A. 
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Let us consider a densely defined self-adjoint operator A. For any g £ D{A), we set 

\\g\\A = [{Ag,Ag) + {g,g)]^/\ (124) 

Then || • \\a is a norm in D{A). As well, we consider {f,g)A — {^1,^9) + if, 9) to be a scalar product specific to 
D{A). 

If /„ € D{A) is a Cauchy sequence in 1 1 • that is, lini„ ^-j-oo ||/m — /nlU = Oj then /„ also is a Cauchy sequence in 
the Hilbert space norm ||/|| = [(/, f)]^^"^- By the completeness of H there is / € "H such that lim„_yoo ||/ — /nil = 0. 
It follows that / necessarily belongs to D{A) (that is, -0(^4) is complete in the || • \ \a norm). Then A = A* is closed 
and denoted A ^ A. (The closcdness is a consequence of D{A) = D{A*) and an assumption that D(A) is dense in "H, 

M-) 

In our text, we have introduced Hamiltonian-type operators H = Hq + V, where Hq is a nonlocal generator of Levy 
noise, via a canonical quantization of the Levy-Khintchine formula. Consequently, we encounter standard quantum 
mechanical domain problems for unbounded operators. This refers as well to H and position-momentum (fairly 
standard !) operators x,p = — iV which form a canonical pair [x,p\ C il, up the presence of h which in various parts 
of the paper has been identified with unity.. 

Let us take Hilbert space vectors f,gEL^ and let (/, g) G C represent a standard scalar product in the complex 
functions space. We need 

if,Hg) = {f,Hog) + {f,Vg) (125) 

to have granted the existence status (c. g. the convergence of involved integrals). In particular, if / is a normalized 
function, we may expect the mean energy value (/, Hf) — (Hq) + (V) here-by defined up to dimensional constants, to 
exist. Since our main objective are domain issues for Hq, let us set V = and pay attention to {f,Hog) exclusively. 

In passing we note that in reference to the semigroup dynamics and the associated jump-type processes, a. well 
developed mathematical formalism of Dirichlet forms associated with noise generators —Hq is available, [isl. ItiI. [t^ . 
The pertinent quadratic forms "live" in a Hilbert space of real- valued functions. This subject matter has been also 
addressed in the physics-oriented literature but exclusively in relation to —A, i. e. in the context of the Wiener noise 
and the emergent semigroup description of the Brownian motion, [2[ . 

We are somewhat guided by basic tenets of that theory, but our considerations will refer directly to the unitary 
(quantum) dynamics and therefore our Hilbert space contains complex valued functions, the real ones forming 
merely a subset in and not a linear subspace, [73{ . 

All symmetric Levy noise generators are Hermitian (actually self-adjoint, {T^I), {f,Hog) = {Hidf,g), on their 
domains od definition. They are also non-negative. Therefore their operator square roots are always well defined as 
self-adjoint operators. 

1/2 

Accordingly, we can pose an auxiliary domain problem for , by considering 

{f,H^g)^{Hl'^f,Hl'^g) (126) 

as a relevant (e.g. non-trivial) part of the induced scalar product (/, p) 1/2 in D{hI'^) C c.f. Eq. (1121. 

The operator domain of interest would consist of all / and g in L^, such that the completion of a corresponding 
linear space in the ||.|| norm would identify a closed dense subspace of L^. We shall not be more elaborate on 

purely mathematical issues and assume such domain existence for granted, for each self-adjoint Hq of Section HI 
separately. 



B. Test model: Quasirelativistic Hamiltonian. 

As a test model let us consider the quasirelativistic Hamiltonian Hm = V — A — m, {h = c = 1). We know that 
solutions of the Salpeter equation 

'^ ^"^tf^ = Wm^ - A - m] V^(x, t) (127) 

have the form ip{'x.,t) = exp(imt) 0(x, t). Therefore, instead of introducing the domain for Hm proper, which is a 
nonnegative operator (with as the bottom generalized eigenvalue), we may pass to an equivalent domain analysis, 
carried out for the strictly positive operator Hm + m — \Jm} — IS., see e.g. our discussion of Section IV. B. 2. 
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This step is advantageous, since if ■(/'(^i^) is a solution of the pseudodiffcrential-Schrodinger equation idtt/j = 
—A + — m]tp, then 0(x, t) = t) exp(— imi) necessarily is a positive energy solution of the free Klein-Gordon 

equation (□ + m2)(/)(x, t) = 0, with □ = -A + d'^/dt^. 

Each scalar positive energy solution 0(x, t) of the free Klein-Gordon equation can be represented in the manifestly 

Lorentz covariant form: 

0(x,t) = J d4fce(-^'^>-")<5(fc^fc''-m2)e(fco)$(fco,k) = ^ J ^$(fco, k)e(-^*'=°+^'^'') (128) 

where k :— {ko, k), fc^fc^ := fcg — k^, $(fc) is a scalar and Q{ko) is the Heaviside function equal to 1 if fcg > and to 
otherwise. Ultimately we are left with fco = Vk^ + 'ti^. 

This representation extends to all solutions of idt(t> ~ V— A -I- m?(j), and upon changing fco — —fco in 0(fco) followed 
by a complex conjugation, to solutions of the time adjoint equation as well. (Side comment: general solutions of those 
pseudodifferential-Schrodinger equations form Lorentz invariant subspaces in the linear space of all solutions to the 
free Klein-Gordon equation.) 

Let us tentatively adopt the following definition of the Klein-Gordon scalar product, [l^, (being independent of 
the specific space-like surface of integration): 

(</>!, </'2)kg i / d''x[fidt^2 -dt'^)cj)2]. (129) 

Given positive energy solutions of the free Klein-Gordon equation (f)i{x),i = 1.2, which we know to solve the 
quasirelativistics equation as well, we realize that (0^ is a Fourier transform of (f>i) 

i(bi,h)KG^ I ^<Pi^2 = (01, V-A + m2</,2) = ((-A + m2)i/40i, (-A + m'Y'^<P2) (130) 
J ko 

We recall that ^/^AT m? is a Hermitian operator in the Hilbert space L^, equipped with a scalar product ((f)i, <f)2) ~ 

J (Px(j)i{x,t)(j)2 ix,t). 

We can now introduce a new positive energy solution for both Klein-Gordon and quasirelativistic Schrodinger 
equations as follows 

^(x) ^ [(-A + m2)i/V](a;). (131) 

At this point we consider as the i^-normalized function in the domain of (m^ — A)^/^, so arriving at the mean 
energy of the quantum system in the state </> (to be compared with Eq. (|126p ): 

((m^ - A)i/40, (m^ - A)i/40) = (</), (m^ - A)!/^^) ^ > 0. (132) 

Accordingly, 

$ = (l//B;)(m2-A)i/40 (133) 

is normahzed, (3>, $) = 1. (We remember that (0, 0) = 1 — >■ (0, 0)kg = -£^0)- 

The normalized function $ may serve as a reference state that gives account of the energy spatial distribution, 
c.f. 0, Hull- The reverse operation obviously reads = ^/E^{m'^ - A)~^/''$. 

Introducing the pdf p(x, t) — |$p(x, t) we have in hands a probability measure. Let 51 be a volume in R^. Then, we 
interpret pe{^) = Jq d'^x |$p as a fraction of the total mean energy that is confined in the volume fl. Compare 
e.g. analogous considerations in the context of the photon wave mechanics in 0, H^l- 

We note in passing that the way we have introduced $ stays in an intimate relationship with the (induced) notion 
of the Newton- Wigncr position operator, [iS, 75] although this notion appears to be irrelevant for our discussion, 
compare also 

In particular, we emphasize that the canonical quantization method adopted by us relies on the standard position- 
momentum pair [xj,pk] C 5ijl, withp = — iV. Accordingly, most familiar uncertainty relations hold true AxAp > 1/2. 
One may merely enter a dispute of what are the states that may sharpen the Heisenberg bound. No new position 
operator notions are here necessary (c.f. past discussions of the covariant position operator notion and the causality 
issue in relativistic quantum theory). 

Our test model discussion allows us to come to a general conclusion. Let us disregard introductory KG equation 
hints (126)-(130). Then, the mean energy-based probability measure definitions (131), (132) can be readily extended 
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(up to dimensional coefficients) to any Levy stable case. To this end one needs to begin with the mean energy 
definition in a given state and replace the quasirelativistic operator \Jm?' — A in defining formulas (131) and (132) by 
a fractional operator |A|''/^ with [i G (0, 2). 

Clearly, any connection with standard (relativistic) wave equations is lost, with a notable exception of the Cauchy 
case /J. = 1, where the D'Alembert link persists. In the latter context, let us assume that the function $(fco, k in Eq. 
(|128p is selected so that imder the to ^ limit, (l/|k|)<i>(|k|, k) has the property to vanish if we let |k| go down to 0. 
Then, all arguments following Eq. (|128p retain their validity, if we r eplace (-A + to^)!/" by |A|i/2" everywhere, for 
n=l,2. 



C. Foldy's synthesis for m > and its m = extension to photon wave mechanics. 

Let us come back to the Salpeter equation in its m > versions. The link of the to > equation with various 
rclatvistic equations has been established long time ago under the name of Foldy's "synthesis of covariant particle 
equations", [7^, see also [t^]. It has been shown that Dirac, Klein-Gordon and Proca equations can be reduced to a 
canonical form 



idtx = \/-^ + m^f3x (134) 

where /3 is a diagonal hermitian matrix. In particular, solutions of the Klein-Gordon equation (□ + m'^)(j)(x,t) = 
yield a two-component wave function x with components x± = (l/\/2)[(— A -I- m'^)~^/^dt4> T «(— A + to^)^/'*(/)]. The 
matrix /3 is identical with the Pauli matrix (73, comprising 1 and —1 on the diagonal. 

The reduction of the Dirac equation to the canonical form is accomplished by the Foldy-Wouthuysen (FW) trans- 
formation. Let ^l> he a solution of 

idtip = {f3m — id ■ ip . (135) 

The canonical form (|134p is obeyed by the four-component x = Uip where U stands for the F-W transformation, 
while /3 comes from the Dirac equation in its traditional block-diagonal form, comprising the 2x2 identity matrix / 
and its negative — /. 

The Proca equation adopted in (76j leads to a canonical form with a 6-component x and /? being being block-diagonal 
with 3x3 blocks / and — /. 

The case of mass to = has been properly addressed quite recently, being associated with attempts to give meaning 
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to the (single) photon wave mechanics, 0,[2j,[25|, see also [57[. The key element there was a reformulation of Maxwell 



equations in terms of the Riemann-Silberstein vector function F(x, t), actually interpreted as the R-S wave function, 
e.g. a solution of the Schrodinger-type equation with a divergence constraint 

idtF = cV X F (136) 
V-F = 0. (137) 

The above dynamical equation can be given more explicit Schrodinger form by invoking spin 1 (3 x 3, no ft involved) 
matrices [Sj, Sk] = iejkiSi, so that c.f. [7| 

idtF =-{S-V)F (138) 

i 

We note that the divergence condition excludes the potentially troublesome fc = mode. Therefore we can safely 
introduce the inverse of the Cauchy operator |A|^/^, while acting upon F such that V • F = 0. 

By introducing the helicity operator A for photons (loosely speaking, a projection of spin on momentum) which is 
a manifestly nonlocal operator, we can rewrite the photon Hamiltonian iif = p • S as follows 

H = c\p\k (139) 

A=^(P-S) (140) 

where p = — iV and wc interpret |p| as |V| — \ ± Aj^/'^. We can pass to the helicity basis (with the corresponding 
eigenvalues ±1), ending with positive energy solutions of the to = variant of the Salpeter equation, in the canonical 
form 



z9x = c|A|i/2/3x 



(141) 
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where /3 is the 6x6 diagonal matrix, encountered before in connection with the Proca equation. Here, x is a 6- 
component vector composed of positive frequency pieces taken form the original R-S vector and its complex adjoint, 
F(+) and F^"^* respectively. Compare a discussion of the helicity /spin spectral issues in Section 2.1 of [3|, see 



also [57 1 . 



In the helicity basis, the photon wave function dynamics is nonlocal and shows all intriguing features (radial expan- 
sion) discussed before in subsection III.F.l. We note that the ex pan ding/contracting behavior has been previously 
associated with the photon energy density F • F* , see e.g. Fig. 3 in [23| . A discussion of involved uncertainty relations, 
from varied perspectives, including the mean energy normalization of the wave function, can be found in 0, |24| . [25j . 
see specifically Eqs. (5.28)-(5.31) of Ref. Q)- Their obvious validity is not that illuminating as a useful indicator 
delocalization in view of rather rapid radial expansion of wave packets. 

By invoking the scalar product considerations of Section 5.1 in we readily arrive at the m = version of the inter- 
play between solutions of the second order wave equations (here, d'Alembert) and those inferred from the Schrddinger 
type equation. Namely, the inversion relation (f> = ^ [rr? — A)~^/"'$ of subsection IV. B has a counterpart in the 
photon wave mechanics. Up to the mean energy normalization, the Landau-Peierls wave function </) = |-A|-i/4$ (c.f. 
section 5.3 in Q) is precisely the outcome of a mapping of solutions of the d'Alembert equation, (^^ — A)F-|- = 0, 
(|128p (set m=0 therein), to I? normalized solutions of Eqs. (138) and /or (141). 

To be more explicit, we note that all previous steps (129) through (130) can be safely repeated if: (i) in view of 
m = 0, we replace (—A -f m^)^/^" by |A|^/^", with n = 1,2, (ii) next extend the interpretation of integration to 
encompass summmation over helicity indices A = 1, —1, and (iii) insert $^(^0, fc) instead of the previous ^(fco, fc)- 

Then the (m | 0) KG scalar product (jl29p assumes the d" Alembertian form and refers to two-component positive 
energy solutons of both d'Alembert and the Cauchy-Schrddinger equation (here we consider the squared norm of ): 
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and the I? normalization of a positive energy solution ((/), 0) = 1 clearly implies (<?!), (f>)d'A — and the normal- 
ization (*, €>) = 1 of 

* = (l/v/^)|A|i/40. (143) 

Now we are exactly at the point where the Landau-Peierls (LP) wave function happens to be introduced. This issue 
was a subject of disagreements in the literature, that ranged form accepting to discard th L-P functions as an object 
of a physical relevance, 0, E^]- On the basis of (142), (143), we know that both $ and (j) arc fairly legitimate objects. 

However if $ is what we interpret to provide a consistent probabilistic interpretation of a single photon quantum 
state (see e.g. also sections 5 and 5 of and [S^]), then we may always reproduce (j) with the squared norm 
{(l),(f>)d'A = Ea/ t^^^*^^ according to (j> = y/E^\A\-^/^<^. We recall that ((^,0) = 1 ^ (0,0)rfM = E^. Hence 
(j)/ \/E^ is d'Alembert-normalized. That can also be accomplished by writing 

-^(0, ^)d'A = (|A|-i/V, \A\-'^U)d'A = 1 . (144) 

We note that often expressed in the literature past prejudices, concerning the LP wave function nonlocality (reaching 
back to W. PauH, (1933), 0), in the light of our discussion must be taken under scrutiny as invalid objections. All 
ingredients of our theory of quantum motion are spatially nonlocal. That is possibly somewhat blurred on the Fourier 
level of description, but fairly obvious while invoking the primary level of of the theory, based on the exploitation of 
non-Gaussian noise (Levy) generators. 



VI. OUTLOOK 



The main message of the present paper is a comprehensive analysis of the intimate relationship between jump-type 
stochastic processes (like e.g. Levy flights) and nonlocal (due to intcgro-differential operators involved) quantum 
dynamics. In the course of the paper we were very explicit in deriving certain formulas (like e.g. those for various 
propagators) which are often presented in the faulty form in the current literature. We have also attempted to clarify 
a number of topics related to so-called fractional quantum mechanics, which as well is not free of drawbacks and 
plainly invalid formulas, [3^ [66|. Such situaton is in part a consequence of serious technical difficulties. Those in fact 
have prohibited the pertinent theory from gaining popularity in the research literature of the 70ties, c.f. Q for earlier 
references. 
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Another serious problem pertains to quite appealing and both customarily and uncritically adopted Fourier rep- 
resentation of the spatially nonlocal quantum dynamics. This is convincingly presented in the literature on various 
aspects of relativistic wave mechanics and that of quantum fields, jH,!!^), where the Fourier analysis has been always 
regarded as a legitimate departure point in the development of quantum theory. 

In the present paper the Fourier representation of analyzed dynamical problems is never a primary, but a secondary, 
derived notion. It does not provide a fully adequate transcription of the original spatially nonlocal formalism, unless 
employed with due care. 

Our input of the theory of random motion, specifically that of jump-type heavy-tailed probability laws, so far 
has not received the status of relevance in the quantum context. Presumably it was the reason of why the inherent 
nonlocality of the quasirelativistic dynamics (take that as an exemplary case) has not been properly addressed. We 
know about one isolated attempt in this direction due to [Toj). 

As discussed in detail, the nonlocal quantum dynamics may lead to strongly enhanced spatial delocalization of 
any "narrow-peaked" initial probability density function. This is not merely a familiar wave packet spreading of 
the standard Schrodinger quantum dynamics. The pertinent delocalization happens to be conspicuously amplified, 
like e.g. in case of simple relativistic physics-motivated models. That is exemplified by a dynamically generated 
multi-modality emerging out of initially unimodal pdfs in the ID case. In the 3D case the quantum dynamics induces 
a radial expansion of any initially localized pdf. The speed of that expansion from the vicinity of the origin, in 
the massless case equals the velocity of light c. This phenomenon makes doubtful a diagnostic utility of standard 
uncertainty relations. 

A point worth mentioning is that of a predominantly lacking "particle" concept of the nonlocal quantum dynamics. 
Doubts pertaining to "what is actually jumping there" have been expressed before Q and refer to the massive Salpeter 
case as well, see also Interestingly, the same "what is jumping" problem persists in the purely classical framework 
of Levy stochastic processes. 

From the jump-processes perspective (which are nonlocal phenomena as well) one should be be very careful while 
interpreting the associated transport phenomena in terms of " particles" . All Levy stable-supported dynamical sce- 
narios are devoid of any genuine "massive particle" connotations. Presently, our view is that it is not a "particle" 
that undergoes a jump- type process, although other opinions have been expressed in the past, [1, Q. 

Since the Salpeter (m > 0) case belongs to our non-Gaussian (and non-Laplacian) inventory, the very "particle" 
position (operator) concept needs to be put under scrutiny, including past vigorous discussions of the causality problem 
in relativistic quantum mechanics and that of of a covariant position operator, or position operator at all in the photon 
context. The "photon position" and related localization concept, appears to be useless in the developed framework 
and turns out to stay in a deep contradiction with the radial expansion of photon wave packets (see below). 

Interestingly, our methods find some support in investigations of the photon wave function (and photon wave 
mechanics). Specifically while combining the time evolution of wave packets with various forms of uncertainty relations 
proposed so far, 0, [Hill, see also j2l| . 
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